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1. INTRODUCTION 

Unit commitment is a significant issue in both the operation of electrical power systems and the 
competitive electricity supply industry and market. It includes two decision-making processes: unit 
scheduling and economic dispatch (ED). The generators scheduling process involves setting the status of 
generating units every hour to be (ON or OFF) of the horizon planning while taking into account the 
constraints of the units' minimum up and down-time. Process of determining the best power generation of 
generation plants (units) to meet power demand and spinning reserve at every hour at the lowest production 
cost within capacity unit limits using economic dispatch goals [1]. 

Because of load variation problem through the day, and for large power systems which contain 
many plants or units, it is not suitable to operate all plants in time. So a method of coordination among these 
plants and generating units must be adopted. The coordination involves in which units will be a start-up and 
the sequence of operating units that must be connected to the network (i.e. ON) and which must be switched 
OFF. Another problem must be solved that it produces electrical power with fulfilling demands in a 
minimum cost [2]. 

Many methods have been used to solve the unit commitment (UC) problem as integer programming 
[3], branch-and-bound methods [4], dynamic programming (DP) [5]-[8], mixed-integer programming [9], 
and lagrangian relaxation methods (LR) [10], [11], and priority list method [12]. Also, there are many 
optimization algorithms suggested to solve unit commitment as whale optimization-differential evolution and 
genetic algorithm (WODEGA) [13], on mixed-integer programming formulations for the unit commitment 
problem [14], solving unit commitment and economic load dispatch problems using modern optimization 
algorithms [15], sested particle swarm optimization (PSO) [16], dynamic formulation of the unit commitment 
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and economic dispatch problems [17]. In this paper the optimal UC solution based on binary grey wolf 
optimizer based on particle swarm optimization (BGWOPSO) for (39 bus IEEE system with 10 units) has 
been presented in discussion and the results compared with related works. 


2. RESEARCH METHOD 

Unit commitment (UC) issues are often used to assess which units must participate in dispatch 
operations at particular times. The number of units, load capacity, startup cost, spinning reserve, and ramp 
rate are all constraints that must be considered. The goal is to create a secure UC system that meets the 
aforementioned constraints at a low total cost. According to their characteristics, power generating systems 
may be listed as base-load, intermediate-load, or peak-load, so the main constraints can be listed as follows: 


2.1. Constraints 

In order to solve UC problem, the constraints must be satisfied which are imposed as operational 
requirements. The constraints can be explained as following [1], [18]: 
a. Power constraint 

The sum of power generation of each generator at time t (Pi) equal to demanded power: 


NJ pt — pt 
Yi=1 Pgi = Po (1) 
where P$ represents active power which demanded at time t and Ng number of generators. 
b. Spinning reserve constraint 

This term is used to characterize the total maximum active generated power from all units in the 
system minus the current load demand. In other words, mathematically we can say: 


A Poi (Max) 2 P$ + PR (2) 


where Peijmax) represents the maximum real power of the unit ( i), and P$ is the reserve of real power at time t. 
c. Power limit constraints 
Pg imin) is is the minimum active power of unit i. 
t 
Pa i(minyS Poa < Pogi(max) (3) 


d. Minimum up time constraint 
Tin is the number of hours that unit (i) is on till time t, and MUT;i is the minimum uptime of unit (i). 


t 
Tion 2 MUT; (4) 
e. Minimum down time constraint 

Tto ff Tepresents off hours of unit (i) till time t, and MDTi is the minimum downtime of unit (i). 


Tfr = MDT, (5) 


In unit commitment, the priority list to make what unit is ON is based on the rate of the cost of full- 
load (y) for a unit which can define as the fuel cost divided on generated power at the maximum value of the 
unit [1]. 


yı = (6) 
gimax 

; = a;Primax + bp t” 7 

Yi l gmax l Pgimax ( ) 


Where a unit that has y; the least with respect to other units will be prioritized to dispatch its power firstly. 

The scheduling of units perhaps does not satisfy the constraint of spinning reserve. Therefore must solve this 

violation that occurred. One approach that found to this purpose is depending on y; as following [1]: 

a) At each time hour, find y; via (7) for each uncommitted unit at hour t, then sort them in an ascending 
form. 

b) Find the requirement of the spinning reserve at t by (2). 
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c) In case the result occurred in point (2) achieves the constraint of spinning reserve, check t, if t < T, go to 
point 1; otherwise, end this process. In the event point 2 is not achieved, one of the uncommitted units 
must be committed that has the least y;from point (1). 


2.2. The economic dispatch 

The main constraint of the operating system is the summation of the active output power of units 
equal to power demand. The input to every unit, shown as fi, where fi refers to the cost rate of operation of 
the unit (i) which is the fuel cost. The output of every unit, Pgi, is the active power generated by the ith unit. 
The total cost represents the summation of the cost of each individual unit which mathematically denoted Fr. 
So the problem of economic dispatch is minimizing Fr. The fuel cost f; of unit (i) at any interval time is a 
function for the power generated of generating unit output Pg; [19]. 


(ff = ci + biPji + a,P3:) (8) 


Where ai, bi, and cj are the fuel cost coefficients. The unit has incremental fuel cost (Ai) which considered the 


slope of the input-output curve and equal an = Ai. So, for the economical partition of load among units 
gi 


within the plant, the norm is that all units must be operating at the same incremental fuel cost. 


2.3. Hybrid grey wolf optimizer base on particle swarm optimization 
2.3.1. Particle swarm optimization 

Particle swarm optimization (PSO) is the algorithm presented by Kennedy and Eberhart [20]. This 
theory was developed relying mainly on the intelligence of the swarm to do this advanced computational 
procedure. In a swarm, every particle is considered a feasible solution to the optimization problem. The 
movement of particles in space depends on the rules of flow. In the case of one of the particles discovering a 
solution, the other particles will approach it. For all individuals, the values of fitness are estimated and that 
will be optimized through the function of fitness. As a PSO is a tool of optimization, it provides an 
exploration process in which residents constantly change their locations over time. Through investigation, 
each particle maintains the best position it found, Prest. In addition, among all the Pres: values gained the better 
fitness is specific by one particle, which is known as the global best, known Gpest [21]. 

The position of the particle and velocity must be updated. Based on the fitness values the personal 
location Phest and global location Gvest are reconstituted as in (9)-(11) [1]: 


Vere = w'Vi cari (Piest E xi) + cari (Ghest T x$) (9) 


xtti = xt + Varn (10) 


i 


Wmax—Wmin 


wt: Inertia weight factor = Wmax — x Iter (11) 


Itermax 


39 66 39 66 


where: “ci, c2: acceleration coefficient”, “x: the location of particle”, “rı, r2: Two numbers taken randomly 
between 0 and 1”. “Vit: The particle’s rapidity (velocity) in i® dimension”. Also where “max and Wmin are the 
initial and final weights and itermax is the maximum iteration number”. 


2.3.2. Grey wolves optimizer (GWO) 

This method relies on grey wolves' behavior in obtaining their food, as grey wolves are at the 
topmost of predators [22]. The hierarchy of grey wolves is the leader in the herd is called alpha (a) and is 
responsible for giving decisions. The second level is beta (£) wolves. The third and lower ranks of the pack 
are delta (6) and omega (w) respectively. Alpha is the fittest solution, both beta and delta will be the second 
and third best solutions respectively [22]. Hunting prey begins first with the process of surrounding prey 
(encircling), which can be modeled mathematically as in (12) and (13) [22]: 


> 


D = |Č. p(k) — Xw(k)| (12) 
Xw(k +1) = Xp(k) — A.D (13) 


where “k: the present iteration”, “A and D: Coefficient vectors”, “X p: The position vector of prey”, 
“Xw: The position vector of the grey wolf”. It can calculate the vectors A and Č as in (14) and (15): 
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A=24.7%-a (14) 
C=2.45 (15) 


e? 


where throughout iterations components “& are decreased from 2 to 0 linearly” and “(77 , 72) are random 
vectors between [0, 1J”. 

Gray wolves can modify their position about the prey [23] and have the capability to distinguish the 
prey's location. The hunting process is usually alpha-guided. To mathematically represent hunting behavior, 
it suggests that alpha, beta, and delta have the best known probable site of the prey [22]. Therefore, the first 
three optimal solutions it was obtained are kept to update their sites in the space proportional with the optimal 
place [24], and the other solutions, including Omega solutions, are neglected [22]. These solutions can be 
determined as in (16) [22]. 


Da = |C-Xa(k)— ŽI, Dg = |C2-Xg(k) — ŠI, Dg = |C. X6(k) — Ž| (16) 


And by the following formulation can calculate the prey position vector with regard to (a, B, and ô) 
respectively [22], [23]: 


Xl = Xa — Ay (Da), Xz = Xg — Az. (Dg), X3 = Xs — As. (Ds) (17) 


where “X7,%2,X3 are a three best solutions (wolves) in the pack at the given iteration k”. By taking the 
average of x1, X2, X3 , it can calculate the best position as the formulation (18) [23]. 


X1 +X2+ X3 
3 


X(k+D= (18) 


As mentioned earlier @ are decreased from 2 to 0 linearly proportional with the number of iterations, the 
algorithm has been updated by the (19) to calculate a [24]: 


2 


“max iter 


G=2-k (19) 


where k: an iteration number, max iter: number of total iterations. 


2.3.3. The hybridization 

Thus, to hybridize the two meta-heuristic methods: GWO and PSO, the exploitation capability in 
PSO and the exploration capability in GWO will be improved to get the strength of the two variables together 
[25]. The first three optimal solutions are modified by the control the exploitation and exploration of the grey 
wolves in the search area by inertia weight factor as in (20) [25]: 


Da = |Ci-Xa(k) —w.X|, Dg = |C-Xg(k) —w.X|, Ds = |C3.X5(k) — w. ŽI (20) 
Then velocity and position equations in (9) and (10) respectively in hybrid GWOPSO will be [25]: 
VE ea! (VÉ + art Ge — xt) + corf (xz — x£) + ard (x3 — x$) (21) 


xtti = xt + pitt (22) 
2.4. Proposed method hybrid binary GWO based on PSO 

Because of the algorithm hybrid GWOPSO is unable to solve the problem of feature selection since 
it is binary problem by its nature, consequently, a binary version of the hybrid PSOGWO must be grown to 
be appropriate for the feature selection problem [24]. Formulation (20) shown the modifying mechanism of 
wolves be a function of three vectors' sites that is x1, x2; x3 and it promotes each wolf to the first three best 
solutions, the updating of position in (18) can be upgraded into (23) [24]: 


. 5 . y, [XL +X + X3 
xkti = f if sigmoid (=a) > rand (23) 
0 otherwise 
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where “ x¥*1: the binary position which updated at kth iteration in dimension d” “rand: a random number 


pulled from the uniform distribution € [1,0]” and sigmoid (a) is denoted as in (24) [24]: 


1. 


sigmoid(a) = EEG] 


(24) 
X1, X2, X3 in (17) are updated and they are calculated by using formulation (25) [24]: 


xt 


f if (x2 + bstep$) 2 1 
0 otherwise 


1 if (xf + bstepł) > 1 
gat I Gp kitaa 25) 
0 otherwise 


xi- : if (x$ + bstep%) > 1 
3 0 otherwise 


x. 


where xå g the position’s vector of (alpha, beta, and delta) wolves in dimension d. 
bstep¢ : a binary step in dimension d that can be denoted as in (26) [24]: 


if cstep gs = rand 


bstep gs = a (26) 


otherwise 


where rand: a random number pulled from the uniform distribution € [1, 0], (d): refers to dimension, 
cstepå gs: Value of the d’s continuous and can be calculated by the next equation [24]: 


1 
ag dpd 
ee CATE p875) 


cstepy gs = (27) 


Thus after update the best three solution positions through (25), in BGWOPSO the exploitation and 
exploration are controlled by an inertia weight factor as in (28) [23]: 


Da = |Ci-Xa (k) —w.X|, Dg = [C2 X p(k) — w.X|, Dg = |C3-X6(k) — w. X| (28) 
Accordingly, the velocity and positions have been updated as in (29) [24]: 

VE = wh (VÉ + crt (x — xi) + Cory (xz — xi) + cará (x3 — xi) (29) 
Note that in (29) the best three solutions x1, x2, x3 are updated according to (25) [24]: 

aftt = ykt1 4 ytti (30) 
where xk*? and V;{*? are calculated based on (23) and (29) respectively [24]. 


2.4.1. Unit commitment program structure using BGWOPSO 
a. Initially, configured particles randomly and input other parameters as iteration number, cl, c2, and w. 

within acceptable limits of generated power (pmin and pmax). 
b. Calculate production cost (fj) and all ED requirements as shown in section (2.3). 
c. Compare the present (f;) for each particle with Pres. If the present (fj) is superior to Pres, then this value is 
Poest otherwise Poest same. Then determine Gres among Prest 
Take the best site in (23) as search agent (individual) current site and modify the site of a, p, 6 wolves 
Take the site of alpha as the final site of the swarm and the alpha result as the best fitness. 
Update the speed of each individual from (29). 
Revise the position of individual x; using (30). 
If the iterations number arrives at the maximum number go to stride (9), else back to stride (3). 
Determine the final cost of all combinations, division of power among the units, and trace the units 
scheduling. 

Engender the best value of Gpest means that it is the best power that can be generated from each unit 

with the lowest total generation cost. Figure 1 presents the above steps through flow chart. 


rood ho e 


' 
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Initialize PSO parameters as: populations number, no. of iterations, c1,c2, W 


o 


No 


Goest (iter) = Gvest (iter-1) Goest (iter) = Gvest 


is Gpest (iter) > Grest (iter-1) 


Take the best site in equation (18) as search agent (individual) current site and modify the site of a,8,5 wolves 


Update particle velocity from (29) and position from (30) 


Display the status Evaluate total cost, 
no. of committed 
units and power 
distribution between 


units 


of committed 
units, load sharing 
between various 
units and total cost 


Figure 1. Flow chart of BGWOPSO algorithm in order to solve UC problem 


3. RESULTS AND DISCUSSION 

The proposed algorithm in this paper was applied on standard test system which is: 39 bus IEEE, 10 
generating units shown in Figure 2. The date of the test system is presented in Tables | and 2 presents the 
daily load for the test system also. The program for this algorithm was built using MATLAB tools 
environment. The optimal parameters of PSO which used to obtain the best result are cj=1, co=1, Wmini=0.4, 
Wmaxi=0.9, and the number of populations is 50 with spinning reserve is 0.1 from load demand. 

The results that arrived of simulating the suggested algorithm based on the test system (IEEE 
39-bus) will be shown in Table 3 (see in appendix) where the units that participate in the operating system 
are recognized in ping color. The final production cost for this case during 24 hours is (563940 $). The 
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results of the IEEE 39-bus test system are compared with related works that using the same standard test 
system under the same conditions (i.e 0.1 spinning reserve and same constraints) as shown in Table 4 (see in 
appendix). It is very clear that the simulation results show the activity of the suggested algorithm BGWOPS 
in the reduction of the operating cost by (0.046939%-0.003294%) respectively, which proves the superiority 
of BGWOPSO over the other optimization approaches. 


39 21 
2 
A 
19 23 
20 36 


9 31 33 2 © 


© 6 CC 


Figure 2. IEEE 39-bus test system diagram 


Table 1. Data of 10 units, 39-bus IEEE system 
Unit No. Prsimax Paimin a($/MW?) b($/MW) c($/h) MUTi MDTi HSCi CSCi CSHi _ISi 


U1 455 150 0.00048 16.19 1000 8 8 4500 9000 5 8 
U2 455 150 0.00031 17.26 970 8 8 5000 10000 5 8 
U3 130 20 0.002 16.6 700 5 5 550 1100 4 -5 
U4 130 20 0.00211 16.5 680 5 5 560 1120 4 -5 
U5 162 25 0.0398 19.7 450 6 6 900 1800 4 -6 
U6 80 20 0.00712 22.26 370 3 3 170 340 2 -3 
U7 85 25 0.00079 27.74 480 3 3 260 520 2 -3 
U8 55 10 0.00413 25.92 660 1 1 30 60 0 -1 
U9 55 10 0.00222 2721 665 1 1 30 60 0 -1 
U10 55 10 0.00173 27.79 670 1 1 30 60 0 -1 
Table 2. Load demand for 10 units system 
Hour 1 2 3 4 5 6 T 8 9 10 11 12 
Demand 700 750 850 950 1000 1100 1150 1200 1300 1400 1450 1500 
Hour 13 14 15 16 17 18 19 20 21 22 23 24 


Demand 1400 1300 1200 1050 1000 1100 1200 1400 1300 1100 900 800 


4. CONCLUSION 

This paper suggested a BGWOPSO algorithm to find solving to the unit commitment. The standard 
test system (IEEE 39-bus test systems) has been used to show the efficiency of the suggested method. From 
the achieved results, it is clear that the BGWOPSO algorithm is succeeded to fulfilling the minimum 
production cost and optimal scheduling for power generation during 24 time hours for the test systems under 
consideration. In terms of computation time and cost, the proposed method (BGWOPSO) outperformed the 
other various technologies in this analysis. The proposed algorithm reduced computing time, allowed fewer 
units to be started up, and had a lower scheduling cost, according to the results. This method can be used by 
power companies to dispatch units, and governments can use the results to promote and implement renewable 
energy power generation. The future work can be used the multiobjective function for minimizing the 
operation cost instead of a single objective function. 
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Table 3. Scheduling of generation and commitment of 10 units system by BGWOPSO 


Table 4. Comparison total production cost between BGWOPSO and other approaches with 


Generation scheduling 


Hour Demand Ul U2 U3 U4 U5 U6 U7 U8 U9 Ul0 
700 455 245 0 0 0 0 0 0 0 0 
750 455 295 0 0 0 0 0 0 0 0 
850 455 370 0 0 25 0 0 0 0 0 
950 455 455 0 0 40 0 0 0 0 0 
1000 455 390 0 130 25 0 0 0 0 0 
1100 455 360 130 130 25 0 0 0 0 0 
1150 455 410 130 130 25 0 0 0 0 0 
1200 455 455 130 130 30 0 0 0 0 0 
1300 455 455 130 130 85 20 25 0 0 0 
1400 455 455 130 130 162 33 25 10 0 0 
1450 455 455 130 130 162 73 25 10 10 0 
1500 455 455 130 130 162 80 25 43 10 10 
1400 455 455 130 130 162 33 25 10 0 0 
1300 455 455 130 130 85 20 25 0 0 0 
1200 455 455 130 130 30 0 0 0 0 0 
1050 455 310 130 130 25 0 0 0 0 0 
1000 455 260 130 130 25 0 0 0 0 0) 
1100 455 360 130 130 25 0 0 0 0 0 
1200 455 455 130 130 30 0 0 0 0 0 
1400 455 455 130 130 162 33 25 10 0 0 
1300 455 455 130 130 85 20 25 0 0 0 
1100 455 455 0 0 145 20 25 0 0 0 
900 455 425 0 0 0 20 0 0 0 0 
800 455 345 0 0 0 0 0 0 0 0 

Units scheduling 
U2 U3 U4 U5 U6 Us U9 U10 


24 


700 

750 

850 

950 

1000 
1100 
1150 
1200 
1300 
1400 
1450 
1500 
1400 
1300 
1200 
1050 
1000 
1100 
1200 
1400 
1300 
1100 
900 

800 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
Hour Demand Ul 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 


E ee E E AE E E EEE E E EE OO E 


The total cost= 563940 $ 


a ee OO 


SODORP RP BRR RRR RR ERR RHR HOO OOO 


© 


ee ee 
SORPR RESP RRR BR BR BeBe RR BR RR OOO 
Orr rr OOOO e FP RFF rr EE A Ba E aA a a 
COFPFRFOCOCOOCORRPH HHH OOSCOSC CO OG 


o 


SCOoOSCOHF COC OC OHHH HCO C OC OCC S 
SecoeoeooCoCoOCoOoC OHH COSCO COCOCOSCSO 
ecoecoeoco oC OC OH COC COC COCSCS 


No. Method Best generation Percentage of production cost 
cost ($) for BGWOPSO to another 

algorithms 

1 New Genetic Approach (NGA) [26] 591715 0.046939 

2 Two Stage Genetic Algorithm (TSGA) [27] 568314.56 0.007694 

3 GA [28] 565866 0.0034036 

4 PSO-LR [29] 565835 0.003349 

5 LR [28] 565828 0.003336 

6 Improved DA-PSO [1] 565807.31 0.003299 

F Discrete Binary PSO [30] 565804 0.003294 

8 Proposed method hybrid BGWOPSO 563940 - 
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